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Abstract. In this paper, the authors first establish the global well-posedness of strong solutions 
of the simplified Ericksen-Leslie model for nonhomogeneous incompressible nematic liquid crystal 
flows in two dimensions if the initial data satisfies some smallness condition. It is worth pointing 
out that the initial density is allowed to contain vacuum states and the initial velocity can be 
i -^h " arbitrarily large. We also present a Serrin's type criterion, depending only on Vd, for the 

breakdown of local strong solutions. As a byproduct, the global strong solutions with large 
initial data are obtained, provided the macroscopic molecular orientation of the liquid crystal 
materials satisfies a natural geometric angle condition (cf. |19j). 
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1 Introduction 

Liquid crystals are substances that exhibit a phase of matter that has properties between those 
of a conventional liquid and those of a solid crystal (cf. [12]). The hydrodynamic theory of 
liquid crystals was first developed by Ericken and Leslie during the period of 1958 through 
1968 (see O [10\ [201 121j). Since then, many remarkable developments have been made from 



both theoretical and applied aspects, however, many physically important and mathematically 
fundamental problems still remain open. In this paper, we consider a simplified Ericken-Leslie 
model for the nonhomogeneous incompressible nematic liquid crystals in two dimensions: 

p t + dw(pu) = 0, (1.1) 

(pu) t + d'w(pu <g>u) + VP = Au - Vd ■ Ad, (1.2) 

divu = 0, (1.3) 

d t + u- Vd = Ad+\Vd\ 2 d, (1.4) 



*This work was partially supported by NNSFC (Grant Nos. 11271306 & 10971171), the Fundamental Research 
Funds for the Central Universities (Grant Nos. 2010121006 & 2012121005), and the Natural Science Foundation 
of Fujian Province of China (Grant No. 2010J05011). 

' Corresponding author (E-mail:jwzhang@xmu. edu.cn). 
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where p : R 2 x [0, oo) — > R + is the density of the fluid, u : M 2 x [0,oo) — > R 2 is the velocity 
field of the fluid, P : R 2 x [0, oo) -> R is the pressure of the fluid, and d : R 2 x [0, oo) -> S 2 
(the unit sphere in R 3 , i.e. |d| = 1) represents the averaged macroscopic/continuum molecular 
orientations. 

Though system (|l.ip — (|1.4p is a simplified version of the Ericksen-Leslie model, but it still 
retains the most interesting mathematical properties without losing the basic nonlinear structure 
of the original Ericksen-Leslie model [HI HU1 EH EI]. Roughly speaking, the system (jl.ip - (jl.4p 
is a system of the nonhomogeneous Navier-Stokes equations for incompressible flows coupled 
with the equation for heat flow of harmonic maps, and thus, its mathematical analysis is full 
of challenges. In particular, if p = Const., then it turns into the following homogeneous system 
which models the incompressible flows of nematic liquid crystal 

u t + u- Vu + VP = pAu - Vd ■ Ad, (1.5) 
divu = 0, (1.6) 
d t + u-Vd = Ad+\Vd\ 2 d (1.7) 

with \d\ = 1. Moreover, if u = in (jl.5[> — (|1.T|) . then it reduces to the following equation for heat 
flow of harmonic maps: 

dt = Ad + \Vd\ 2 d, \d\ = l. (1.8) 

There has been a lot of literature on the mathematical studies of (|1.5j) - (|1.7j) and (11.8p . see, 
for example, [HI H3 [23l [Ml EH [261 E31 [35] and [3 El [H El [32], respectively. In the following, 
we briefly recall some related mathematical results of the liquid crystal flows. In a series of 
papers, Lin [23] and Lin-Liu [24\ [25] initiated the mathematical analysis of (|1.5p — (|1.7p in 1990s. 
More precisely, to relax the nolinear constraint \d\ = 1, they proposed an approximate model 
of Ericksen-Leslie system with variable length by Ginzburg-Landau functionals, that is, the 
equation (jl.7p with \d\ = 1 is replaced by 

d t + u-X7d = Ad+i I (l-\d\ 2 )d. (1.9) 

In [23l [23] , the authors proved the global existence of classical and weak solutions of (jl.5p , (|1.6p , 
(jl.9p in dimensions two and three, respectively. The partial regularity of suitable weak solutions 
was also studied in [25]. However, as pointed out in |24j . the vanishing limit of e — > is an open 
and challenging problem. Indeed, in contrast with (11. 9p . it is much more difficult to deal with 
the nonlinear term | Vd\ 2 d with \d\ = 1 appearing on the right-hand side of (jl.4p or (jl.7p from the 
mathematical point of view. In two independent papers |13j and [26] . Hong and Lin-Lin- Wang 
showed the global existence of weak solutions of (]1.5p - (]1.7p in dimensions two, and proved that 
the solutions are smooth away from at most finitely many singular times which is analogous to 
that for the heat flows of harmonic maps (see [31 E2] ) . The global existence of smooth solution 
with small initial data of f|1.5[) (jl.7|) was also proved [261 S3] an d [351 122] in dimensions two and 
three, respectively. 

For the approximate nonhomogeneous equations (jl.ip - (jl.3p and (|1.9p . the global existence 
of weak solutions with generally large initial data was proved in [281116]. and the global regularity 
of the solution with strictly positive density was studied in [8]. As aforementioned, the nonlinear 
term \Vd\ 2 d with \d\ = 1 will cause serious difficulty in the mathematical analysis of liquid 
crystal flows. Recently, Wen and Ding [3l] established the local existence and uniqueness of 
strong solutions of (|l.ip -( fT~4j) in the case that the initial density may contain vacuum states 
(i.e. po > 0). Moreover, if the initial density has a positive lower bound (i.e. po > p > 0) which 
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indicates that there is absent of vacuum initially, the global strong solutions with small initial 
data was also obtained in [34] . 

As that for the density-dependent Navier-Stokes equations (see [H[27]), the possible presence 
of vacuum is one of the major difficulties when the problems of global existence, uniqueness and 
regularity of solutions are involved. Therefore, in the present paper we aim to investigate the 
global regularity of (jl.4j) when the initial density may contain vacuum. 

We consider the Cauchy problem of (ll.ip - (ll.4p with the following initial data: 

(p,u,d)(x,0) = (p ,u ,d )(x) for x G K 2 , (1.10) 

and the far-field behavior at infinity: 

(p, u, d)(x, t) — )■ (p, 0, e) as \x\ — > oo, t > 0, (1-H) 

where p > is a given positive constant and e G § 2 is a given unit vector (i.e. |e| = 1). 

To state our main results, we first introduce the definition of strong solutions of (|l.ip - (|1.4p . 
(fLTUjl and (fTTTj) . 

Definition 1.1 A pair of functions (p,u,P,d) is called a strong solution of (jl.ip - (|1.4p . (|1.10p 
and (fTTTJ) on R 2 x [0,T], if p(x,t) > for all (x,t) el 2 x [0,T], 

' p-pGC([0,T];H 2 (R 2 )), p t G L°°(0, T; H^R 2 )) 

u G C([0, T];H 2 (R 2 )) n L 2 (0, T; H 3 {R 2 )), 
< y/put £ L°°(0,T;L 2 (IR 2 )), u t G L 2 (0, T; H\R 2 )), (1.12) 

VP G C([0, T]; L 2 (M 2 )) n L 2 (0, T; H 1 ^ 2 )), 
k Vd G C([0, T]; ff 2 (M 2 )), d t G L°°(0, T; i7 1 (R 2 )) n i 2 (0, T; H 2 (R 2 )), 

and (p,u,P,d) satisfies (TTTD -dT^ll a.e. on R 2 x (0,T]. 

Then, our first result concerning the global strong solutions with small data can be stated 
in the following theorem. 

Theorem 1.1 Assume that the initial data (po,UQ,do) satisfies 

Po > 0, (po - P, uo, Vd) G H 2 (R 2 ), divn = 0, |d | = 1, 
Auo - VP - Vd • Ad Q = pl /2 g for some (VP ,g) G L 2 (R 2 ). 

Then for any given < T < oo, there exists a unique global strong solution (p, u, P, d) of 
(fj~T|) -([L4l. (fTTTO]) and (fLlTI) on R 2 x [0,T], provided 

exp (2 (|| P y\ |l! 2 + ||V(io|| 2 2 )) ||Vd || 2 2 < 1. (1.14) 

It is worth mentioning that the smallness condition ()1.14|) stated in Theorem 11.11 implies 
that (po, uo) can be arbitrarily large if || V^oIIl 2 i s chosen to be suitably small. This is analogous 
to the one in [33]. Moreover, as a result, we see that the strong solution to the Cauchy problem 
of nonhomogeneous Navier-Stokes equations (i.e. d = Const.) with large initial data, which 
may contain vacuum, exists globally on R 2 x [0, T] for all < T < 00. Thus, Theorem 11.11 also 
generalizes the result due to Huang- Wang [14"] . 



3 



The proof of Theorem 11.11 is mainly based on a critical Sobolev inequality of logarithmic 
type which was recently proved by Huang- Wang (cf. [15]) and is originally due to Brezis-Wainger 
[2] (see also [291 EI]). However, it is remarkable that the arguments in [15] actually depend on 
the size of the domain considered and cannot be applied directly to the case of the whole space. 
Thus, some new ideas have to be developed. The main difference lies in the proof of Lemma 
13.31 where, instead of llp 1 ' 2 ^*!!^ and 1 1 /o 1 / 2 -^ • Vn|| i 2, we use the material derivative ||p 1//2 n||/ / 2 
for some technical reasons. We also note here that the strictly positive far-field condition p > 
plays an important role in our analysis. The strongly nonlinear terms |V(i| 2 (i and Vd • Ad in 
(|1.2p and (|1 .4[) will also cause some additional difficulties. 

For the generally large initial data, it is still an interesting and open problem whether the 
strong solution blows up or not in finite time. In [26] and |14j . the authors proved respectively 
that the following blowup criteria for the two-dimensional equations of (|1.5p - (|1.7|) : 

f' T f T 

li™ / (Nit* + II Vd||i*) dt = oo and lim / \\Vd\\ L ~dt = oo, (1.15) 
Jo Jo 

where < T* < oo is the maximal time of the existence of a strong solution to f)1.5[) ()1.7|) . 
Motivated by the proofs of Theorem II. 1[ we can prove the following mechanism for possible 
breakdown of strong solutions, which is a natural extension of the ones in |26l 114] . 

Theorem 1.2 Assume that < T* < oo is the maximal time of the existence of a strong 
solution to (|l.ip ~( [L~4"j) . (jl.lOp and (jl.lip with generally large initial data (po,uo,do) satisfying 
(fTT3j) . Then, 

lim / \\Vd\\' L rdt = oo (1.16) 
1 Jo 

for any (r, s) satisfying 

111 . 

- + -<-, 2<r<oo. (1.17) 

r s 2 

Theorem 11.21 implies that for any < T < oo if the left-hand side of (|1.16|) is finite, then 
the strong solution of ([TI|> - ([r4]> . (fTTUJ) and (fTTTTJ) will exist globally on M 2 x (0,T). 

Based on a frequency localization argument combined with the concentration-compactness 
approach, Lei-Li-Zhang [19] recently proved the following interesting rigidity theorem for the 
approximate harmonic maps. 

Proposition 1.1 ( 119\, Theroem 1.5]) For given positive constants < Cq < oo and < e < 1, 
assume that d : M 2 — > S 2 satisfying Vd G H 1 (M. 2 ) with ||Vd||^2 < Co and d% > e. Then there 
exists a positive constant 5q £ (0, 1), which depends only on Co and e, such that 

||W||£ 4 <(1-<5 )||V 2 < 2 , (1.18) 

which particularly implies 

|| Ad + |Vd| 2 d|| 2 2 > | (|| Ad||£ a + || Vd\\%) . (1.19) 

As an immediate consequence of Theorem 11.21 and Proposition II. 1[ we can remove the 
smallness restriction (|1.14p on the initial data and prove the following existence theorem of 
global strong solutions with large initial data, provided the macroscopic molecular orientation 
of the liquid crystal materials satisfies a natural geometric angle condition. This extends the 
Lei-Li-Zhang's result (cf. [19J) to the case of nonhomogeneous incompressible liquid crystal flows 
with initial vacuum. 
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Theorem 1.3 Let e% = (0,0,1) G S 2 and let do3 be the third component of do. Besides the 
condition (11 . 13|) in Theorem assume further that 

d 03 > e and d - e 3 G L 2 (M 2 ) (1.20) 

/io/ds /or some uniform positive constant e > 0. T/ien /or any < T < oo, there exists a unique 
global strong solution (p,u,P,d) of (fTTj) - flUD ■ (fTTTU]) and (Oil on R 2 x [0,T]. 

The rest of the paper is organized as follows. In Sect. 2, we state some known inequalities 
and facts which will be used later. The proof of Theorem 1 1 . 1 1 will be done in Sect. 3, based on 
the local existence theorem and the global a priori estimates. In Sect. 4, we outline the proof 
of Theorems 11.21 and 11.31 



2 Preliminaries 

In this section, we list some useful lemmas which will be frequently used in the next sections. We 
first recall the well-known Ladyzhenskaya and Sobolev inequalities (see, for example, [171 ^1)- 

Lemma 2.1 For f G H 1 ^ 2 ), it holds for any 2 < p < oo that 

|| 4 < V^||/|M|V/|| i2 , (2.1) 
\ LP <C(p)\\f\\% p \\Vf\\ 1 L - 2/p , (2.2) 

where C(p) is a positive constant depending on p. In addition, if f e W 1 *^ 2 ) n H 2 (M. 2 ) with 
p > 2, then there exists a universal positive constant C such that 

II/IIloo < C\\f\\ m , P < C\\f\\ m . (2.3) 

We will also use the following Poincare type inequality, which shows that the velocity u 
actually belongs to L 2 -space even that the vacuum states may appear. 

Lemma 2.2 Let p > be a given positive constants. Assume that g — p G L 2 (M 2 ) n -L°°(R 2 ) 
with g(x) > 0, Vv G L 2 (R 2 ) and ^fgv G L 2 (R 2 ). Then, 

\\v\\ L 2 < C{p, \\g - p\\&nL°°) (\\p 1/2v \\l* + UVullia) , (2.4) 

where C(p, \\g — pW^HL 00 ) i> s a positive constant depending only on p, \\q — p\\l 2 an d \\q~ p\\l°°- 
Proof. Indeed, by virtue of Holder and (|2.2j) . we have for any q > 2 that 

\v\ 2 dx = jM 2 dx-j\g-p)\ V \ 2 dx 

< CWg^vWl, +C( \g- p\idx) ( / \v\ 2 ^^dx\ 
^ru~ ii ~n s / ii 1/2 m2 , ii i|2(g-l)/gnv7 i|2/g" 

<c{p, \\q - p\\ L 2nL°°) [\\s f nU + pw ^ " q \\Vv\\ '? 



< C(p, \\g- p\\ L 2riL°°) (\\Q 1/2 v\\h + ll^llla) + |lMliaj 
which proves (|2,4p immediately. □ 

Next, to improve the regularity of the velocity, we need to use the following estimates of 
the Stokes equations (see, for example, [TT| [18]). 
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Lemma 2.3 Consider the following stationary Stokes equations: 

-AU + VP = f, div?7 = in R 2 . 

Then for any f £ W m,p {R 2 ) with m E Z + and p > 1, there exists a positive constant C, 
depending only on m and p, such that 

\\V 2 U\\ w ^p + ||VP||w/m,P < C\\f\\w».p- (2-5) 

To estimate the L 2 -norm of the gradient of the velocity, we shall apply a critical Sobolev 
inequality of logarithmic type which was prove by Huang- Wang (cf. [15]) and is originally due 
to Brezis-Wainger [2] (see also [29\ [3Tj ) . This is the key tool for the proofs of Theorems I1.1H1.31 

Lemma 2.4 For q > 2 and < s < t < oo, assume that f G L 2 (s, t; H l (R 2 ))nL 2 (s, t; W 1,9 (M 2 )). 
Then there exists a positive constant C(q), independent of s,t, such that 




In the case that the lower bound of the density is nonnegative, the local existence of strong 
solutions to (|l.ip - (jl.4p . (jl.lOp and (jl.lip was proved in [34|. Indeed, in [M] the authors only 
considered the case of smooth bounded domains, however, as pointed out in [7], the similar 
procedure also works for the whole space by means of the standard domain expansion technique. 
For simplicity, we quote the following local existence theorem of strong solutions without proofs. 

Lemma 2.5 Assume that the conditions of Theorem hold. Then there exists a positive time 
< To < oo such that the Cauchy problem (jl.ip - (jl.4p . (jl.lOp and (jl.lip admits a unique strong 
solution on R 2 x (0,T ). 

3 Proof of Theorem 11.11 

Assume that the conditions of Theorem 11.11 hold . Let < T* < oo be the first blowup time of a 
strong solution (p, u, P, d) to the Cauchy problem lll.lj) - (|1.4|) . (ll.lOp and (|1.11|) . In order to prove 
Theorem 11.11 it suffices to prove there actually exists a generic positive constant < M < oo, 
depending only on the initial data (po,uo,do) and T*, such that 

£(T)= sup ( ||p - p\\ H 2 + \\u\\ H 2 + ||Vd||#2 + \\p 1/2 u t \\ 2 L 2 + \\d t \\ 2 H i ) 
o<t<T v 7 

+ [ T (IMIffa + \\Vd\\ 2 H3 + \\u t \\ 2 Hl + lldtl^a) dt<M (3.1) 
Jo 

holds for any < T < T*. So, by the local existence theorem (see Lemma I2.5P it can be easily 
shown that the strong solution can be extended beyond T*, which gives a contradiction of T*. 
Hence, the strong solution exists globally on M? x [0, T] for any < T < oo. The proof of 
Theorem 11.11 is therefore complete. 

The proof of (|3.ip is based on a series of lemmas. Throughout the remainder of the paper, 
for simplicity we denote by C a generic constant which depends only on the initial data and T*, 
and may change from line to line. 
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First, it is easy to see from the method of characteristics and (jl.ip that for every < T < T* , 

< p(x,t) < \\p \\ L °° for all (x, t) € M 2 x [0, T). (3.2) 

Moreover, multiplying (jl.ip by q\p — p\ q ~ 2 {p — p) with q > 2, integrating it by parts over (0, t), 
and using the divergence- free condition f j 1 . 3 j) , we find that 

II (P-P)(*) II w = ||po-/5||w for Vte[0,T]. (3.3) 
In view of (|l.ip - (|1.4p . we have the following standard energy estimates. 
Lemma 3.1 For every < T < T* , one has 

sup f (\p l/2 u\ 2 + |V<i| 2 ) dx + 2 j I (\Vu\ 2 + \Ad+\Vd\ 2 d\ 2 )dxdt 

< J (\p 1/2 u\ 2 + |Vd| 2 ) (x,0)dx 4 E . (3.4) 

Proof. Multiplying (I1.2P by u in L 2 and integrating by parts, by (I1.3P we know that 

^— y p|u| 2 (ix + J \Vu\ 2 dx = - J (u ■ Vd ■ Ad) dx. (3.5) 

Due to the fact that \d\ = 1, multiplying (jl.4p by (Ad+\Vd\ 2 d) in L 2 , we obtain after integrating 
the resulting equations by parts over M 2 that 

-T- I \^d\ 2 dx+ [ \Ad+\Vd\ 2 d\ 2 dx 
2dt J Jq 

= J (u-Vd- Ad)dx + J (\Vd\ 2 d-d t + \Vd\ 2 u-Vd-d) dx 

= J (u- Vd- Ad)dx+ l - J (\Vd\ 2 d t \d\ 2 + \Vd\ 2 u- V\d\ 2 ) dx 

= j (u ■ Vd • Ad) dx, (3.6) 

which, combined with ()3.5|) . immediately leads to (|3.4p . □ 

To be continued, we need the following key estimates on || V 2 (i|| i 2( 0i T;L 2 )- 
Lemma 3.2 Assume that the initial data satisfies 

ex P (2(||p; / \o||i2 + ||Vdo||i 2 ))||Vd || 2 o2 < ^, (3.7) 
then it holds for every T G (0, T*) that 

rp 

sup ||Vd||| 2 + / \\V 2 d\\ 2 L2 dt < — . (3.8) 
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Proof. After integrating by parts, we easily deduce from the identity \d\ = 1 that 

J \Ad + \Vd\ 2 d\ 2 dx = J (\Ad\ 2 + |V(i| 4 ) dx - 2 j \Vd\ 2 (d ■ Ad)dx 

= J (\Ad\ 2 -\Vd\ 4 )dx. (3.9) 

On the other hand, integration by parts, together with the divergence-free condition (jl.3p . gives 

J (u ■ Vd ■ Ad) dx = - J {d j u i d i d k d j d k + i);/ i)/'^ dx 

= ~ j (d j u i d i d k d j d k ^ dx < ||V«|| L2 ||Vd|| 2 4, (3.10) 

where and in what follows the repeated indices denotes the summation over the indices. 
Putting (|3.9p . (I3.10P into (|3.6j) and recalling the fact that 

\\Ad\\ 2 L2 = \\V 2 d\\ 2 L2 , 

we obtain 

^IIV* + HV 2 d|ll 2 < \\Vu\\ L 4Vd\\% + ||Vd||| 4 , 



which, combined with (J2JJ) and the Cauchy-Schwarz inequality, yields 

lj t \\™\\h + HV 2 d|| 2 L2 < V^||Vn|| i2 ||Vd|| L2 ||V 2 d|| L2 + 2|| W||i 2 ||V 2 d||i 2 



< (2||Vd|| 2 2 + ||V 2 d|| 2 2 + 2||Vu|| 2 2 ||Vc/||| 2 . (3.11) 

It follows from ^7} that 

||Vdo||i 2 <e 2 ^||Vd |li 2 <^, 

and thus, by the local existence theorem and the continuity argument we see that there exists 
a T\ > such that for any t G [0, Ti], 

l|Vd||| 2 < ±. (3.12) 

Set 

f = sup{T | (I3TT2D holds}. 
Then it follows from (l3~TT]) ~ (l3~T2l that for any t G [0,T), 

|||Vd|| 2 2 + ||V 2 d|| 2 2 <4||Vn|| 2 2 ||Vd|| 2 2 , 
which, together with Gronwall's inequality and (|3.4p . leads to 

l|Vd||£a+jf \\V 2 d\\ 2 L2 dT <exp (±J ||Vu||| 2 dr^ ||Vd ||| 2 



<e 2Eo ||V(io|li2 <^. (3.13) 

Combining (|3.4p . (13. 13ft with the continuity argument immediately implies that f)3. 12|) holds for 
all < f < T*, and thus, the proof of <^M) is finished. □ 

By Lemmas 12.41 and 13.21 we can now derive the estimates of ||Vti||^2 and ||V 2 d||£ 2 which is 
the most important step among the proofs. 

Lemma 3.3 For every < T < T* , one has 

sup (Hlri + llvdl&i + Hllla) 

0<t<T 

+ £ (||V 2 u|| 2 2 + \\p 1/2 u\\ 2 L2 + \\dtfjjx + ||Vd||^ 2 ) dt < C, (3.14) 
which particularly gives 

rT 

\p 1/2 u t \\ 2 L 2 < C. (3.15) 



/ 

Jo 
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Proof. Let / = ft + u ■ Vf denote the material derivative. Also set 

1 c)d c)d 

M(d) 4 Vd^Vd- -|Vd| 2 I 2 , (Vd® Vd)ij 4 _ - . _ , l<i,j<2, 

& CJJL 'i CJJL j 

then it is easily seen that 

Vd- Ad = div(M(d)). 

To prove (13.14|) . multiplying (11. 2p by ut and integrating it by parts over M 2 , we deduce 



1 d 
2di' 



VU||^2 + \\P 1/2 U\\ 2 L 2 



= — J divM(d) • Ufdx + J pu ■ Vu ■ udx 

= — J M(d) : Vudx - J M(d) t : Vudx + J pu-Vu- udx 

<jj M{d) : Vudx + \\\p X / 2 u\\l. + \\\Vd t \\ 2 L2 + C(\\uf Laa + ||Vd||£ 00 )||V«||£ 3 , 

where we have also used (|3.2H and Cauchy-Schwarz inequality. As a result, 

j t \NAh + \\p 1/2 n\\ 2 L 2 

< 2— / M(d) : Vudx + -||Vd t || 2 2 + C(\\uf Loa + ||Vd|||oo)||Vu|| 2 2 . (3.16) 

Next, one easily obtains from (|1 .4j) that 

|l|Vd|| 2 2 + (K||| 2 + ||V 2 d||i 2 ) < C y* (M 2 |Vd| 2 + |Vd| 4 ) dx 

<C(hlli- + l|Vd|| 2 00 )||Vd|| 2 , 2 . (3.17) 
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To deal with the term ||Vc^|| 2 2 ° n the right-hand side of (|3.16p . we first apply V to both 
sides of (jl.4j) to get that 



Vd t - VAd = -V(u • Vd) + V(|Vd| 2 d), (3.18) 

from which it follows that 

j t \\V 2 d\\ 2 L2 + (\\Vd t \\ L2 + \\VAd\\l 2 ) 



< J (|V(|Vd| 2 d)| 2 + |V(n- Vd)\ 2 )dx 

<C J (\Vd\ 6 + \Vd\ 2 \V 2 d\ 2 + \Vu\ 2 \Vd\ 2 + \u\ 2 \V 2 d\ 2 ) dx 

< C {Ml- + llVdlli-) (|| Vu||£ a + ||V 2 d|| 2 2 ) , (3.19) 
where we have used (|2.ip and (|3.4p to get that 

||v< 6 < ciivdii^iivdiil.iiv^n 2 , < ciiwiilooliv 2 ^. 

Using (|2.ip and (|3.4p again, we have 



M(d) : Vudx < \\\Vu\\ 2 L2 + C\\Vd\\ 2 L2 \\V 2 d\\ 2 L2 < \\\Vu\\ 2 L2 + C^V 2 d\\ 2 L2 . 

Taking this into account, multiplying (|3.19l) by 2C± + 1, and adding the resulting inequality, 
(|3,16p and (|3.17p together, we obtain after integrating the resulting inequality over (s,t) with 
< s < t < T that 



(||V< 2 + ||W|| 2 , 1 )(*) + J [\\ p V^l 2 + \\ dt f Hl + \\^ d f H ^ dT 

< C (||V«|| 2 2 + HVdH^) (s) + c£ (\\uf Loo + HVdllloc) (||Vu|| 2 2 + \\Vd\\ 2 Hl ) dr, 
and consequently, 

(||Vu||| a + ||Vd|| 2 Hl ) (t) + f Up^uWl. + \\d t \\ 2 m + \\V 2 d\\ 2 Hl ) dr 



<C{\\Vu\\ 2 L2 + ||Vd||^x) (s)exp (cj {\\u\\ 2 Loo + llVdHicc) dr^j . (3.20) 
Clearly, it remains to estimate ||(it, Vd)||i,°o. To this end, let 

Ht)±e+ sup (||V< 2 + ||V< 1 )(t)+ r(||// 2 «||| 2 + ||dt|| 2 , 1 + ||V 2 d|| 2 rl )dr. 

0<T<t Jo J 

First, in view of ([Ml) and (|3T^|) - (|57^|) . we have 

\\u\\ 2 L2 < C (\\p 1/2 u\\ 2 L2 + ||Vu|| 2 2 ) < C (1 + ||Vn||i 2 ) . (3.21) 

Next, using Lemma [2T3l (|2.ip . (|3.2p and (|3.4p . we deduce from Holder and Cauchy-Schwarz 
inequalities that 

||V 2 n|| L2 <C(|H| i2 + ||Vd-Ad|| L2 ) 
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< C (||p 1/2 n|| L2 + ||V 2 d||| 2 + ||V 3 d|| L2 ) , (3.22) 
which, combined with (|3.8p . yields 

f l|V 2 n||i 2 dr < C jT (||p 1/2 n||i 2 + ||V 2 d||£ 2 + ||V 3 d||i 2 ) dr 

< C sup ||V 2 d||| 2 + C f (\\ P l/2 u\\ 2 L 2 + ||V 3 d||| 2 ) dr. (3.23) 

Thus, recalling the definition of <3?(T) and using f|3.4[) . (|3.2ip and (|3.23p . we infer from 
LemmaEl]that for any < s < t < T < T*, 

W u \\h(s,t;L°°) < C (l + \\u\\ 2 L2 {Sjt . H i) In (e + \\u\\ L 2 {Stt . w i,4))) 
<c(l+ ||Vu|| 2 2(M;L2) In (e + ||n|| L2(M;H2) )) 

< C (l + ||Vu|| 2 2(M;i2) In (e + ||V 2 n|| i2(Sit;i2) )) 

< C (l + ||Vn|| 2 2(S)t;i2) In (C<&(t))) . (3.24) 

In a similar manner, by (|2.6p and (|3.4p one has 

HVd|| 2 2(Sit;ioo) < C (l + ||Vd|| 2 2(M;H1) In (e + ||W|| L2(Sjt;WM) )) 

< C (l + ||V 2 d||i 2(M;L2) In (e + ||V 2 d|| L2(M;Hl) )) 

<c(l + ||V 2 d|| 2 2(M;i2) ln(C$(t))). (3.25) 

For any < s < t < T < T* , putting (I3T241 and ([335]) into (JOty gives 

<&(*) < C$( S )exp{c 2 (||V U ||| 2(M;L2) + ||V 2 d||| 2(M;L2) ) ln(d$(t))} 

< CT*(.) [C 1 $(t)] C2 ( l|Vu|l ' 2 (^ 2 ) +l|V2d|l - 2 (^- 2 )) . (3.26) 
It follows from (|3.4p and f|3.8[) that there exists a positive constant 5 > such that 

C 2 (||Vn||| 2(T _, jT;L2) + ||V 2 d||| 2(T 
which, inserted into (|3.26j) . leads to 



1 



$(T) < C$(T - 5) [Ci$(T)] 1/2 < -$(T) + C7$ 2 (T - 5), 

so that 

$(T) < C(T)$ 2 (T - 5). (3.27) 

As a result of (|3.27p . we see that $(T) is bounded for any < T < T* since the local 
existence theorem indicates <J?(T — 5) < oo for any < T < T*. This, together with (|3.2ip and 
U332J), finishes the proof of (f3"7i~4"D . 
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Furthermore, recalling the definition of material derivative (i.e. " "), one gets from (|2.ip . 
P3D and (EHIjl that 

/ \\p 1/2 ut\\ 2 L 2dt< [ \\p 1/2 u\\ 2 L2 dt+ [ ( p\u\ 2 \Vu\ 2 dxdt 
Jo Jo Jo J 

<C + C I \\u\\ 2 L 4Vu\\ 2 Li dt 
Jo 

r T 

<C + C \\u\\ L 2\\Vu\\ 2 L 2\\V 2 u\\ L 2dt 
Jo 

<C + C [ \\V 2 u\\ 2 L2 dt < C, 
Jo 

which immediately proves f|3. 15j) . The proof of Lemma 13.31 is therefore complete. □ 
Next, we proceed to estimate ||/3 1//2 nt|| i 2 and ||Vdt[|j^a. 

Lemma 3.4 For every < T < T* , one has 

0<t<T v " J Jo 



sup (ll/^Vll^ + IKH^) + f (\\Vu t \\ 2 L2 + \\d tt \\ 2 L2 + \\V 2 d t \\ 2 L2 )dt<C, (3.28) 

<t<T v J Jo 

and moreover, 

sup (\\ufg2 + \\Vd\\ 2 H2 ) + f 1 {\\Vu\\ 2 wlA + \\V 2 d\\ 2 H2 ) dt < C. (3.29) 

<t<T JO 



-T 

0<t<T ' Jo 

Proof. Differentiating (|1.2|) with respect to t gives 

putt + pu ■ Vu t — Aut = —pt(ut + u ■ Vu) — put ■ Vu — VPt — divM t , 
which, multiplied by ut in L 2 and integrated by parts over M 2 , results in 



2~dt 



--"P^lli. + liv^nl, 



= -/«(«. + «■ V«) ■ - / ^ • V« • t*fe + J M t : Vu t d X 

±h + h + h (3.30) 

2 

where ^4 : B = ajjfoy for A = (0^)2x2 and -B = (6^)2x2- 
».i=i 

We are now in a position of estimating the right-hand side of (|3.30p term by term. First, 
using (jl.ip and integrating by parts, by Lemma 12.11 (|3.2p and (|3.14j) we deduce 

I\ = J (/m • V|«(| 2 + pu ■ V(u ■ Vu ■ ut)) dx 

< C j (p|w||^||Vut| + p\u\\Vu\ 2 \u t \ + p\u\ 2 \V 2 u\\u t \ + p\u\ 2 \Vu\\Vu t \) dx 

< C ^||u||Loo||p 1/2 u i || L 2||Vn t || i 2 + UnllLoollVtill^llp 1 / 2 ^!!^ 
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+C (j|u||ioo \\p 1/>2 Ut\\ L 2 1| V 2 u\\ L * + H|~ II Vnii^nv^n^) 

<^||V« t ||i 2 +C(l + ||pV2 Ut ||4 2 + || V 2 u| |2 2 ^ 

where we have used Cauchy-Schwarz inequality and the following estimate due to Lemma 12.11 
and (^Tlj) : 

\\u\\ L °° < C\\u\\ w i,4 < C (\\u\\ m + ||Vn||^ 2 ||V 2 u||^ 2 2 ) < C (l + ||V 2 n||^ 2 2 ) . (3.31) 
Due to ([Ml), ([22]) and ([33]), we have 

Klli»<Cr([|p 1 / a « t [|^ + [|V« t ||^) > (3.32) 
and thus, by (|2.1|) . (|3.2|) and (|3.14|) the second term I2 can be bounded as follows: 

\h\ < C J p\u t \ 2 \Vu\dx 

<C\\Vu\\ L 4u t \\ L 4p l l 2 u t \\ L 2 
^CllVull^llV^ll^llnJ^llVn.ll^llp 1 ^!!^ 

< C (HV^lUa + ||p 1/2 n t || L2 ) \\v 2 u\\% 2 \\p 1/2 Mv 

< \\\s?iH\\b + c { 1 + W^Mh + IIp 1/2 ^II1 2 ) • 

Finally, it is easily seen from (|2.3p and (|3. 14[) that 

\h\<C\\Vd\\ L ~\\Vd t \\ L2 \\Vu t \\ L 2 

<\\\Vu t \\l 2 +C(l + \\V 3 d\\l 2 )\\Vd t \\l 2 . 

Substituting the estimates of Ii,l2 and I3 into (|3.30p . one obtains 

j t \\p 1/2 ut\\h + \\Vu t \\h 

< C (l + \\V 2 u\\h + Wp^utfv) + C (1 + \\Vdf H2 ) \\Vdt\\l 2 . (3.33) 
To estimate ||Vdt||£,2, we differentiate (|1.4|) with respect to t to get 

dtt- Adt = (\Vd\ 2 d-u- Vd) t , 
and hence, using Lemma I2TT1 and (|3.14p . we deduce after direct calculations that 

^IIVdtHi, + (Rt|| 2 2 + ||Ad t ||| 2 ) 

(|Vd| 2 |Vd t | 2 + |Vd| 4 |d t | 2 + |^| 2 |Vd| 2 + |u| 2 |V(i t | 2 ) dx 

< c (nvdiiloo + iMiioo) ||vd t || 2 2 + c[| vd[|i« ||dt||| 2 + c\\vd\\l4u t \\l A 

< C [\\Vdf H2 + \\u\\ 2 H2 ) ||V^||i 2 + d (||W||^ 2 + ||p 1/2 ^||i 2 + ||Vu t || 2 2 ) , (3.34) 
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where we have also used (|3.32[) and the following estimate due to (|2.ip . (|2.3p and f|3. 14|) : 
||Vd[|i«, < C (||Vd|| 4 L4 + ||V 2 d|| 4 L4 ) < C (1 + ||V 2 ci|| 2 2 ||V 3 d||| 2 ) < C (1 + ||V 3 d||| 2 ) . 
Now, multiplying ()3.33|) by 2C\ + 1 and adding it to ()3.34p . we see that 

j t (\\p 1,2 ut\\h + ||V^||| 2 ) + (||V« t ||| 2 + ||AcZ t ||| 2 + \\d tt \\ 2 L2 ) 
< C (1 + HVdllk + ||u||^,) (1 + ||V^|| 2 2 ) + CWp^utWl^ 

which, combined with (|3.14p . (|3. 15|) and Gronwall's inequality, leads to ()3.28j) . since the com- 
patibility condition stated in (|1.13|) 9 implies that (p l / 2 u t ){x, 0) € L 2 (R 2 ) is well defined. 
Using d2U, <£T2§, (I3~T41) and (I3T28]) . we have by ([23]) that 

||V 2 n|| i2 < C ( 1 1 yoxtt 1 1 z, 2 + \\p u ' Vu||x,2 + \\Vd ■ A<i|| L2 ) 

< c (\\ P l i 2 u t \\ L 2 + IM| L 4||Vu|| L 4 + livdii^nv 2 ^!^) 

< C (l + ||V 2 n||^ 2 2 + IIV^H^ 2 ) . (3.35) 
Similarly, one also infers from (|3.18p that 

||v 3 d|| L2 < c {\\vd t \\ L 2 + ||v(« • vd)|| L2 + ||v(|Vd| 2 d)lM 

< c (i + \W\\ L 4v 2 d\\ L , + \\Vu\\ L 4Vd\\ L A + ||vd||i« + livdn^nv 2 ^!^) 

^cfi + iiv^ii^ + liv^n 172 



\ L 2 T || V U|| i2 

from which, ()3.35p and Young's inequality, we arrive at 



sup (||u|| H 2 + \\Vd\\ H 2) < C. (3.36) 

0<t<T 

Using (I2.5P and Lemma ET] again, we deduce from (|3.2p . (|3.28j) and ()3.36p that 

/ ||V 2 n|| 2 : 4 < C I {\\pu t \\ 2 L i + \\pu ■ Vu\\ 2 L 4 + \\Vd- Ad|| 2 4 ) dt 
Jo Jo 

<C J (\\ P 1/2 u t \\h + ||Vn t ||| 2 + \\u\\ A H2 + ||W||^ 2 ) dt 

<C + C [ \\Vut\\ 2 L 2dt < C. (3.37) 
Jo 

Moreover, by virtue Lemma l2.lt (I3.28P and (I3.29P we infer from (ll.4p that 

f T \\v 2 d\\ 2 H2 dt < c [ T (\\d t \\ 2 H2 + \\ u ■ vd\\ 2 H2 + \\\vd\ 2 d\\ 2 H2 ) dt 

Jo Jo 

<C [ {\\d t \\ 2 H2 + hll^UVdll^ + UVdllfl* + ||Vd||f, 2 ) dt 
Jo 

<C, 

where we have used the following Moser's type calculus inequality (see |30J ) that for f,g G 
H S (R 2 ) with s > 2, 

\\fg\\H' < C (\\f\\ L °o \\g\\ H s + < C||/||h.[|^||h.. 

This, together with (I3.36P and (|3.37p . leads to (I3.29f> immediately. □ 
The last step is to estimate the first and second order derivatives of the density. 
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Lemma 3.5 For every < T < T* , one has 

sup (||Vp|| H i + \\p t \\m) + [ \Wu\\ 2 m dt < C. (3.38) 

0<t<T Jo 



Proof. Differentiating (jl.ip with respect to X( (i = 1,2), multiplying the resulting equation by 
|Vp| p ~ 2 <9i/9 with p > 2, and integrating it by parts over R 2 , we obtain after summing up that 

|||V/C» < c\\vu\\ LX \\v P \\i P < c\\vu\\ wl ,4Vp\f LP , 

which, combined with f|3.29j) and Gronwall's inequality, yields 

\\Vp\\ p LP <C\\Vp \\ p LP exp(c \\Vu\\ w iAdt\ < C, Vp>2. (3.39) 
Similarly, by (|3.38|) we also deduce from (jl.ip that 

j t \\v 2 p\\h < c\\Vu\\ L ~\\v 2 p\\i 2 + c\\v 2 u\\ L 4v P \\ L 4v 2 p\\ L 2 

< C7||Vu|| W i,4 (l + ||V 2 p||i 2 ) , 

so that 

||V 2 p||| 2 <C(l+ HVV0III2) exp [C J \\Vu\\ w iAdt) < C. 

As a result of (pHHD . (pTBUD and ([OOJl . one easily gets from (JTT]| that ||pt||#i < C. 
Finally, it follows from (123]) . (l3~2jl . (|3Ti4|) . (EOHjl and (137381) that 

||Vn||^2 < C(||/)tif ||^i + ||pu • Vu||#i + || Vd • Ad||#i) 
< C (||p 1/2 ^|| i2 + ||Vp|U4||u t |U4 + ||Vu t || L2 ) 

+c {\H 2 H2 + \\u\\ L ~\\Vp\\ L 4Vu\\ L 4 + \\Vd\\ 2 H2 ) 

<C7(l + ||Vu t || L3 ), 
where we have also used Lemmas 12.11 and 12.21 Consequently, 



(3.40) 



f \\Vu\\ 2 H2 dt<C + C [ \\Vu t \\ 2 L2 dt < C, 
Jo Jo 



which, together with ([3739]) and (|3T5D|) . proves ([3738]) . □ 

Collecting all the estimates in (j3.2|) . (|3.3p and Lemmas 13. 1H3.5I together . we arrive at (|3.ip . 
and hence, the proof of Theorem 11.11 is complete. 

4 Proofs of Theorems 11.21 and 11.31 

This section is concerned with the proofs of Theorems 11.21 and 11.31 We first prove Theorem 11.21 
by using contradiction arguments. So, to do this, we assume otherwise that 

rT 



T limJ \\Vd\\ s Lr dt<M Q <oc (4.1) 
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with any (r, s) satisfying (|1.17p . 

We begin the proof with the observation from the proof of Theorem 11.11 that, to remove 
the smallness condition ()1.14p and to obtain a global strong solution with generally large initial 
data, it suffices to achieve the estimate of || V 2 d||^2(o,T;L 2 ) f° r an Y < T < T*. Moreover, it 
follows from d3T9]) and (J33J) that for any < T < T*, ' 

T 







V 2 d\\ 2 L2 dt= [ \\Ad+ \Vd\ 2 d\\ 2 L2 dt+ [ \\Vdf Li dt 
Jo Jo 



1 



T 



~2 Eo + J H V *^- ( 42 ) 

Therefore, to bound ||V 2 d|| i 2( 0iT;i 2), we only need to deal with ||V(i| \l a (o,t-,l a )- This will be 
done in the following. 

On one hand, assume that (r, s) satisfies 

- + - < - with 4 < r < oo. (4.3) 
r s ~ 2 ~ ~ y 1 

Then, using Holder inequality, (|2.ip and (|3.4p . we find 

||Vd||l 4 = ||Vd||| 4 ||Vd||| 4 

< c\\ vd\\%-^ r ~ 2) || vd\0 r ~ 2) || Vd|| L2 1| V 2 d|| L2 
<C||Vci||2( (, - 2) ||V 2 d|| L2 

<l\\V 2 d\\ 2 L2 + C(l + \\Vd\\i r ). (4.4) 
Thus, putting (|4.4p into (|4.2p and using (|4.ip . we obtain 

[ T ||V 2 d||| 2( ii < E + C f T (1 + ||Vd||ir) dt < C, (4.5) 
Jo Jo 

provided (r, s) satisfies (|4.3p . 

On the other hand, assume that (r, s) satisfies 

- + - < - with 2 < r < 4. (4.6) 
r s 2 

Then, by virtue of (|2.2p . (|3.4p and Holder inequality we find that 



||Vd|| L4 < C||V<7 a ||Vd|| 



LP 



<C||Vd||i7 a ||Vd||^ /p ||V 2 d||^- 2)/p 

< C\\Vd\\l7^\V 2 d\\fr 2)/p (4.7) 



with 

2r (4 — r)p 

<p<oo, a = -j- ^ G 



r - 2 r ' 4(p - r) V ' 2 

As a result of (|4.7p . we have by Young's inequality that 

livdiil^ciivdii^iiv 2 ^! 4 ^- 2 ^ 
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2 

since direct calculations give 



< I||v 2 d||| 2 + c\\vd\\f- a)p/{p - 2a{p - 2)) 

<il|V 2 d||| 2 + C(l + ||Vd||10, (4-8) 



4(1 — a)p _ 2r < 



p - 2a(p - 2) r - 2 

Thus, putting f|^78|) into P~2j) , by (|4TT|) we also obtain (|4"3|) . provided (r, s) satisfies (|4T6|) . 

To conclude, we have proved that there exists a positive constant C, depending on the 
initial data, T* and Mq, such that for any < T < T*, 



[ \\V 2 d\\ 2 L2 dt < C, 
Jo 



provided (14. lj) holds. With the help of this and (I3.2p - (|3.4I) . following the arguments in the proofs 
of Lemmas l3.3H3.5l we arrive at (|3.ip . which, combined with the local existence theorem (see 
Lemma 12 . 5(> . implies the solutions can be extended beyond T*. This immediately leads to a 
contradiction of T*, and hence, the proof of Theorem 11.21 is complete. □ 

Proof of Theorem \1.SX In fact, by applying the maximum principle to the equation of d$ 
(i.e. the third component of d), we have 

inf d 3 (x,t) > inf d 03 > e, V t > 0. 

xGM 2 xgK 2 

So, it follows from (|3.4p and (|1.19p that for any < T < oo, 

[ T (||Ad||| 2 + ||Vd||£ 4 ) * < C(E ,£) C \\Ad + \Vd\ 2 d\\ 2 L2 dt < C{E ,e). (4.9) 
Jo Jo 

This, combined with Theorem 1 1 . 2 1 with r = s = 4, implies that the strong solutions of (|l.ip - (|1.4p . 
(jl.lOP and (jl.lip exist for all T > 0. The proof of Theorem 11.31 is thus finished. □ 
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